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Tables and charts are given of mathematical functions related to the 
'potential of a line of point charges. The use of these functions is illustrated 
by applications to semiconductor resistivity measurements and to calculations 
of the base resistance of point-contact transistors. 

INTRODUCTION 

The "method of images" is a simple but elegant technique for the 
solution of problems in potential theory. Formal solutions to problems 
involving plane boundaries and point sources can usually be written 
down immediately from obvious symmetry considerations. The solutions 
thus obtained, however, are often infinite series whose convergence is 
unsatisfactory for numerical calculations. 

Many techniques have been devised for improving upon the con- 
vergence of these infinite sums. It has generally been assumed that the 
mathematician or engineer must select one of these techniques and apply 
it to his own special problem. The object of this article is to make such 
special treatment unnecessary in most cases. To this end, the potentials 
of certain simple image systems are tabulated. 

The use of the tabulations is illustrated by numerical solutions of 
some problems in semiconductor engineering. Transistors and varistors 
have become famous because of their non-ohmic properties. But solu- 
tions to ohmie flow problems can be useful as first approximations in the 
design of semiconductor devices and can be very accurate approxima- 
tions for suitably arranged measurements of the equilibrium resistivity 
of semiconductors. 1 '" Plane boundaries occur-as the physical boundaries 

105 



106 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1955 

of the semiconductor and as junctions between regions of different con- 
ductivities and conductivity types. Point sources are the simplest 
idealizations of point contacts. 

The four-point method to be described for the measurement of semi- 
conductor resistivity is used in geophysics for measuring the resistivity 
of the earth. The mathematical content of this article is applicable to 
problems other than ohmic flow. The differential equations for ohmic 
flow are identical with those of electrostatics, heat conduction, and the 
hydrodynamics of an ideal fluid. In one of the examples given below, 
the Coulomb energy of an ionic crystal is calculated. 

SUPERPOSITION OF AN INFINITE NUMBER OF POINT SOURCES 

The strength q of a point source will be defined so that the potential 
due to the source is q/r, where r is the distance from the source. Sources 
of negative strength are often called "sinks", but this special designation 
will not be used in the discussion to follow. All the cases of ohmic flow to 
be considered here involve point sources of current at a plane surface 
and it will be shown that 

«-i5 (l) 

where / is the current into the plane surface of material of resistivity p. 

According to the principle of superposition, the potential at a given 
point due to some configuration of fixed sources is the sum of the po- 
tentials that would be established at that point by the various sources 
making up the configuration. The computations for many problems 
could be accomplished easily if one could have a table of the potential 
of a line of equally-spaced point sources of equal value q. But a divergent 
expression is obtained when one attempts to calculate the potential of 
such an arrangement of sources by superposing the q/r potentials of 
the sources. One way to avoid this divergence is to abandon the attempt 
(implicit in the choice of q/r as the potential of a source) to make the 
potential at infinity equal to zero. Potentials with respect to some ar- 
bitrary point would be finite and could be tabulated. 

However, having the potential at infinity equal to zero simplifies the 
superposition of source systems. All of the problems that could be solved 
with the potential of a line of point sources can be solved with either of 
two related arrangements that permit setting the potential at infinity 
equal to zero. Tabulations of the potentials of these arrangements will 
be given. 
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The first scheme that will be considered consists entirely of point 
sources and is illustrated in Fig. 1. It consists of a line of equally spaced 
point sources of equal strength q, and a parallel line of opposite sources, 
except that one of the positive sources and the corresponding negative 
source have been omitted. The potential at the point P (where the 
omitted positive source would have been placed) is to be considered and 
may be written 

V(P) = q - M{\) (1) 

a 

where M is a dimensionless function of the ratio X = d/a of the distance 
d between the lines of sources to the spacing a of the charges in the lines. 
It is obvious that M (X) must vanish at X = 0. 

The second scheme is illustrated in Fig. 2. It consists of a line of 
equally spaced point sources embedded in a line source of equal and 
opposite strength per unit length. The potential at the point Q, lying on 
a perpendicular to the line of charges, erected at one of the point sources 
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Fig. 1 — Arrangement of sources Fig. 2 — Source distribution corre- 

corresponding to the function M. The sponding to the function A'. The po- 

potential at point P is (q/a) M(d/a). tential at point P is (q/a) N(d/a). 
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Table I — 


Values of the 


Functions M(X) and N(X) 


X 


•1/(X) 


if to 


X 


if(X) 


A'(X) 


0.01 


0.000120 


90.557677 


0.52 


0.278275 


0.105086 


0.02 


0.000481 


41.943612 


0.53 


0.287536 


0.097637 


0.03 


0.001081 


26.087274 


0.54 


0.296876 


0.090740 


0.04 


0.001921 


18.328465 


0.55 


0.306292 


0.084353 


0.05 


0.003000 


13.773673 


0.56 


0.315780 


0.078434 


0.06 


0.004317 


10.803665 


0.57 


0.325338 


0.072947 


0.07 


0.005872 


8.729459 


0.58 


0.334961 


0.067860 


0.08 


0.007662 


7.209018 


0.59 


0.344647 


0.063140 


0.09 


0.009686 


6.053671 


0.60 


0.354392 


0.058761 


0.10 


0.011943 


5.151023 














0.61 


0.364193 


0.054696 


0.11 


0.014432 


4.430064 


0.62 


0.374047 


0.050921 


0.12 


0.017150 


3.843793 


0.63 


0.383952 


0.047416 


0.13 


0.020096 


3.359908 


0.64 


0.393903 


0.044159 


0.14 


0.023266 


2.955502 


0.65 


0.403900 


0.041133 


0.15 


0.026660 


2.613904 


0.66 


0.413937 


0.038320 


0.16 


0.030273 


2.322701 


0.67 


0.424014 


0.035705 


0.17 


0.034105 


2.072472 


0.68 


0.434127 


0.033273 


0.18 


0.038151 


1.855945 


0.69 


0.444273 


0.031012 


0.1(1 


0.042410 


1.667423 


0.70 


0.454451 


0.028907 


0.20 


0.046877 


1.502384 














0.71 


0.464658 


0.026949 


0.21 


0.051550 


1.357196 


0.72 


0.474890 


0.025127 


0.22 


0.056426 


1.228910 


0.73 


0.485148 


0.023431 


0.23 


0.061501 


1.115110 


0.74 


0.495426 


0.021852 


0.24 


0.066773 


1.013798 


0.75 


0.505725 


0.020381 


0.25 


0.072237 


0.923312 


0.76 


0.516041 


0.019012 


0.26 


0.077889 


0.842254 


0.77 


0.526373 


0.017736 


0.27 


0.0S3726 


0.769448 


0.78 


0.536718 


0.016548 


0.28 


0.089746 


0.703889 


0.79 


0.547074 


0.015440 


0.29 


0.095942 


0.644722 


0.80 


0.557441 


0.014409 


0.30 


0.102312 


0.591213 














0.81 


0.567816 


0.013447 


0.31 


0.108852 


0.542726 


0.82 


0.578196 


0.012551 


0.32 


0.115557 


0.498711 


0.83 


0.538582 


0.011715 


0.33 


0.122425 


0.458690 


0.84 


0.598970 


0.010936 


0.34 


0.129450 


0.422244 


0.85 


0.609360 


0.010210 


0.35 


0.136629 


0.389006 


0.86 


0.619750 


0.009532 


0.36 


0.143958 


0.358654 


0.87 


0.630138 


0.008900 


0.37 


0.151432 


0.330903 


0.88 


0.640523 


0.008311 


0.38 


0.159048 


0.305500 


0.89 


0.650904 


0.007761 


0.39 


0.166801 


0.282222 


0.90 


0.661279 


0.007248 


0.40 


0.174687 


0.260868 














0.91 


0.671647 


0.006770 


0.41 


0.182703 


0.241262 


0.92 


0.682007 


0.006323 


0.42 


0.190844 


0.223244 


0.93 


0.692358 


0.005906 


0.43 


0.199107 


0.206671 


0.94 


0.702698 


0.005518 


0.44 


0.207486 


0.191417 


0.95 


0.713027 


0.005155 


0.45 


0.215979 


0.177364 


0.96 


0.723344 


0.004816 


0.46 


0.224582 


0.164411 


0.97 


0.733647 


0.004499 


0.47 


0.233289 


0.152462 


O.UX 


0.743935 


0.004204 


0.48 


0.242098 


0.141434 


0.99 


0.754209 


0.003928 


0.49 


0.251005 


0.131248 


1.00 


0.764466 


0.003671 


0.50 


0.260006 


0.121836 














1.01 


0.774706 


0.003431 


0.51 


0.269098 


0.113135 


1.02 


0.784928 


0.003206 
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Table I — Values of the Functions M(X) and N(X) — (Continued) 



X 


■l/(X) 


A'(X) 


X 


M(\) 


A'(X) 


1.03 


0.795132 


0.002996 


1.54 


1.280952 


0.000100 


1.04 


0.805316 


0.002801 


1.55 


1.289715 


0.000093 


1.05 


0.815480 


0.002618 


1.56 


1.298447 


0.000087 


1.06 


0.825624 


0.002447 


1.57 


1.307149 


0.000082 


1.07 


0.835746 


0.002287 


1.58 


1.315821 


0.000077 


1.08 


0.845847 


0. 002 138 


1.59 


1.324464 


0.000072 


1.09 


0.855924 


0.001999 


1.60 


1.333077 


0.000067 


1.10 


0.865979 


0.001869 














1.61 


1.341660 


0.000063 


1.11 


0.876010 


0.001748 


1.62 


1.350214 
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1.12 


0.8S601S 


0.001634 


1.63 


1.358739 
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1.13 


0.896000 


0.001528 


1.64 


1.367234 
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1.14 


0.905958 


0.001428 


1.65 


1.375700 


0.00004S 


1.15 


0.915S90 


0.001336 


1.66 


1.384137 


0.000045 


1.16 


0.925797 


0.001249 


1.67 


1.392544 


0.000042 


1.17 


0.935677 


0.001168 


1.68 


1.400923 


0.000040 


1.18 


0.945531 


0.001092 


1.69 


1.409273 


0.000037 


1.19 


0.955358 


0.001022 


1.70 


1.417594 


0.000034 


1.20 


0.96515S 


0.000956 














1.71 


1.425886 


0.000032 


1.21 
1.22 
1.23 
1.24 
1.25 
1.26 
1.27 
1.28 
1.29 


0.974930 
0.984675 
0.994392 

i.oohinii 

1.013740 
1.023371 
1.032974 
1.042547 
1.052091 


0.000894 
0.000836 
0.000782 
0.000731 
0.000684 
0.000640 
0.000599 
0.000560 
0.000524 


1.72 
1.73 
1.74 
1.75 
1.76 
1.77 
1.78 
1.79 
1.80 


1.434150 
1.442385 
1.450593 
1.458772 
1.466923 
1.475046 
1.483141 
1.491208 
1.499248 


0.000030 
0.000028 
0.000026 
0.000025 
0.000024 
0.000022 
0.000020 
0.000019 
0.000018 


1.30 


1.061606 


0.000490 


1.81 


1.507260 


0.000017 








1.82 


1.515244 


0.000016 


1.31 


1.071091 


0.000459 


1.83 


1.523202 


0.000015 


1.32 


1.0805-17 


0.000429 


1.84 


1.531132 


0.000014 


1.33 


1.089973 


0.000401 


1.85 


1.539036 


0.000013 


1.34 


1.099369 


0.000376 


1.86 


1.546912 


0.000012 


1.35 


1.108735 


0.000351 


1.87 


1.554762 


0.000011 


1.36 


1.118072 


0.000329 


1.88 


1.562585 


0.000010 


1.37 


1.127378 


0.000308 


1.89 


1.570381 


0.000010 


1.38 


1.136654 


0.000288 


1.90 


1.578151 


0.000009 


1.39 


1 . 145900 


0.000270 








1.40 


1.155115 


0.000252 


1.91 


1.585895 


0.000009 








1.92 


1.593612 


o.ooooos 


1.11 


1.164300 


0.000236 


1.93 


1.601304 


o.ooooos 


1.42 


1.173455 


0.000221 


1.94 


1.608970 


0.000007 


1.43 


1.182580 


0.000207 


1.95 


1.616609 


0.000007 
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1.191674 


0.000193 
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1.624224 
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1.20073S 


0.000181 
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1.631812 


0.000006 
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1.209772 
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1.639376 
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1.218775 
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1.646913 
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0. 000(H) 1 


1.52 
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0.000114 


2.40 
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0.000000 


1.53 


1.272161 
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2.50 


2 000718 


0.000000 
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may be written 



V(Q) = 2 N(\) 

(X 



(2) 



if we continue to denote by a the spacing between the point sources and 
now take d to be the distance between the point Q and the line. As X 
becomes infinite, N goes to zero. 

The potential of a line of point sources along a perpendicular to the 
line erected at one of the point sources can be taken at 



V 



-iri- 

a La 



MCk) 



(3) 



where X = d/a and d is the distance from the line of sources. 

Numerical values of M and N can be obtained with sufficient accuracy 
for most purposes from the charts in Figs. 3 and 4. More accurate values 
of M and N are given in Table I for various values of X. 

A useful relation exists between the functions M and N. It is 

M (X) + N{\) = \ + 2 In X - 2 In 2 4- 2C, (4) 

X 

where C is Eider's constant, 0.577215665. That is 

M (X) + N(\) = I + 2 In X - 0.231863031 
X 

The methods of computing M and N are described in the appendix. 



SEMICONDUCTOR 




Fig. 5 — Linear four-point probe. 
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FOUR-POINT RESISTIVITY MEASUREMENTS 

Resistivity measurements are regularly being made on germanium with 
a four-point probe. The basic assumptions are that the resistivity is 
uniform and that the potential distribution can be accurately approxi- 
mated by that calculated for ohmic flow. The results of calculations for 
some geometries have been presented by L. B. Yaldes, along with a 
discussion of the conditions for the validity of the assumption of ohmic 
flow. 1 In brief, it is required that modulation of the conductivity by 
injected minority carriers be kept negligible. It is also necessary that 
the probes used to measure potential difference be sufficiently conducting 
to pass the current drawn by the voltage measuring instrument. Both of 
these requirements can usually be met, in the case of germanium, by 
abrading the surface or by subjecting the contact between probe and 
germanium to an electric discharge, i.e., "forming" the contact. 

The prototype geometry for four-point resistivity measurements is 
given in Fig. 5, which is taken from Reference 1. The potentials at 
probes 2 and 3 are 

V 2 = $ - q 



St s 2 + s 3 

(5) 



Vz = 



S t + S-, S 3 



where q is the strength of a source corresponding to the current /. By 
By considering a hemispherical surface of infinitesimal radius r centered 
on one of the current probes, one sees that 

/ = - X area X E = - 2^ %, 
P p r- 

or 

q = Ip/2ir (6) 

where p is the resistivity and E = q/r" is the magnitude of the electric 
field. 
The potential difference AV between probes 2 and 3 is 

\$i <S 3 Sz -\- S3 Si + S2/ 

It will be assumed in the following examples that Si = S2 = S3 = s 
in which case 

AV = q/s (8) 
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so that 

P = 2va(AV/I). (9) 

When the four-point probe is applied to a solid that is not approxi- 
mately semi-infinite, a potential difference will still be observed when a 
current flows, so we will define an "apparent resistivity" po by 

po - 2ts(AV/I) (10) 

The true resistivity will be given by 

P = po/C.D. (11) 

where CD. is the correction divisor for the particular case. 

One of the simplest and most useful departures from the semi-infinite 
geometry is the slab or slice of finite thickness, free of any conducting 
coating on its faces (and resting on a non-conducting support during 
the measurement). Fig. G shows the position of the probes and the 
configuration of sources required for calculating the correction divisor. 
The required configuration can be obtained by superimposing two of the 
arrangements shown in Fig. 1 (one with opposite sign) i.e., upon the 
source and sink used for the semi-infinite geometry. The additional 
potential at one of the probes due to the additional sources is 

£- [M(2s/2w) - M(s/2w)] 

and the additional potential at the other probe is equal and opposite, 
so that the additional potential difference is 

q - \M(2s/2w) - M(s/2w)] 
w 

In comparison to the potential difference for the semi-infinite case, 
we find that AV is "too high" by a factor 



CD. -fff+4 

s w 



» &)-<&!) 



= l + s - 

w 



"';-"•£ 



)1 



(12) 



This expression is true in general, but is most convenient for slices that 
are thick in comparison to the spacing of the probes, since the correction 
terms are then small and errors in reading M from Fig. 3 do not matter 
much. The correction divisor may be given in terms of the function A r 
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through use of equation (4). The result, 



CD. = - 

w 



2\n2 + N[±)-N(- 

\2w/ \w 
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(13) 



is most useful for thin slices. In fact, if the slice thickness w is less than 
half the probe spacing s, CD. « (2 In 2) s/w and equations (10) and (11) 
give 

p « ttwAV/I In 2 £3 4.53wAF/7 (for w « Y 2 s) (14) 

The correction divisor is given in a table in Fig. (> for a number of ratios 
of probe spacing to slice thickness. For larger values of this ratio than arc 
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Fig. 6 — Linear probe on infinite slice with nonconducting faces (or on edge of 
semi-infinite slice). 
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found in the table, equation (14) may be used with sufficient accuracy 
for all practical resistivity measurements. 

In general, a large correction divisor is desirable because it indicates 
that a relatively large potential difference is to be measured. For this 
reason, the use of the four-point probe on thin slices should be and, in 
fact, is found in practice to be quite satisfactory. However, if the slice 
is provided with a conducting coating on one side, as in intermediate 
stages of point-contact transistor fabrication, the measurement is apt 
to be inaccurate. The correction divisor for this case is 



CD. = 1 + 



S -\ 2 M 
w [_ 



2m 



- M 



w) \4w 



(15) 



or 



CD. = 



IV 



ns) + "v5l- 3A " 



8 

2w 



(1(0 



Fig. 7 illustrates this configuration and contains a table of the correction 
divisor. It will be noted that the correction divisor is quite small if the 



V 



< 



i 



\i 





^j^*~— - 


*f 




*Z~^*~- 


CONDUCTING 
" BOTTOM 




s/w 


CD. 


s/w 


CD. 


s/w 


CD. s/w 


CD. 


0.1 

0.2 


.9993 

.9948 


0.5 
0.8 


.9329 
.7960 


1.0 

1.5 


.6833 2.0 
.4159 5.0 


.2283 
.0034 



Fig. 7 — Linear probe on infinite slice with conducting bottom face. 
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Corrcctic 


n Divisor for Conducting Case 










L/s 








0.1 


0.2 


0.5 


1.0 


2.0 


5.0 


10.0 


0.0 


0.034 


0.124 


0.481 


0.811 


0.962 


0.9971 


0.9996 


0.1 


0.03 


0.124 


0.48 


0.81 


0.96 


0.997 


1.0001 


8 °- 2 


0.03 


0.125 


0.48 


0.81 


0.96 


1.002 


1.0064 


- 0.5 


0.04 


0.125 


0.49 


0.83 


1.01 


1.08 


1.090 


w 1.0 


0.04 


0.142 


0.56 


1.03 


1.34 


1.48 


1.497 


2.0 


0.066 


0.22 


0.95 


1.84 


2.46 


2.72 


2.765 


5.0 


0.146 


0.55 


2.35 


4.58 


6.12 


6.78 


6.894 



Correction Divisor for Non Conducting Case 
L/s 








0.1 


0.2 


0.5 


1.0 


2.0 


5.0 


10.0 


0.0 


2.000 


1.9661 


1.8764 


1.5198 


1.1890 


1.0379 


1.0029 


1.0001 


0.1 


2.002 


1.97 


l.ss 


1.52 


1.19 


1.040 


1.004 


1.0017 


0.2 


2.016 


1.98 


1.89 


1.53 


1.20 


1.052 


1.014 


1.0094 


s 0.5 


2.188 


2.15 


2.06 


1.70 


1.35 


1.176 


1.109 


1.0977 


w 1.0 


3.009 


2.97 


2.87 


2.45 


1.98 


1.667 


1.534 


1.512 


2.0 


5.560 


5.49 


5.34 


4.61 


3.72 


3.104 


2.838 


2.795 


5.0 


13.863 


13.72 


13.32 


11.51 


9.28 


7.744 


7.078 


6.969 


10.0 


27.726 


27.43 


26.71 


23.03 


18.56 


15.49 


14.156 


13.938 



Fig. 8 — Semi-infinite slice — linear probe parallel to edge. 

spacing between probe points is much larger than the slice thickness. 
The two-point measurement to be described below seems preferable to 
the four-point measurement in the case of thin slices with a conducting 
coating, because less effort is required to maintain an accurate small 
spacing between two points than between four points. 
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The two cases that have been treated thus far were also given by 
Valdes, but have been mentioned here as simple examples of the use of 
the M and N functions. An example of a more complicated problem 
is that of the linear four-point probe on a semi-infinite slice, parallel 
to the edge of the slice, as shown in Fig. 8. It is useful to note that the 







Correction Divisor for Nonconducting Edge 







s/w 





0.1 


0.2 


o.s 
1 . 1333 


1.0 


2.0 


5.0 


10.0 


« 


0.0 


1.4500 


1.3330 


1.2555 


1.0595 


1.0194 


1.0028 


1.0005 


1.0000 


0.1 


1.4501 


1.3331 


1.2556 


1.1335 


1.0597 


1.0198 


1 .0035 


1.0015 


1.0009 


0.2 


1.4519 


1.3352 


1.2579 


1 . 1364 


1.0637 


1.0255 


1.0107 


1.0084 


1.0070 


0.5 


1.5285 


1.4163 


1.3476 


1.2307 


1 . 1648 


1 . 1263 


1 . 1029 


1.0967 


1.0939 


1.0 


2.0335 


1.9255 


1.8526 


1.7294 


1.6380 


1.5690 


1.5225 


1.5102 


1.5045 


2.0 


3.7236 


3.5660 


3.4486 


3.2262 


3.0470 


2.9090 


2.8160 


2.7913 


2.7799 


5.0 


9.2815 


8.8943 


8.6025 


8.0472 


7.5991 


7.2542 


7.0216 


6.9600 


6.9315 


10.0 


18.5630 


17.7886 


17.2050 


16.0944 


15.1983 


14.5083 


14.0431 


13.9199 


13.8629 





Correction Divisor for Conducting Edge 







s/w 





0.1 


0.2 


0.5 

0.8667 


1.0 

0.9405 


2.0 


5.0 


10.0 


. 


0.0 


0.5500 


0.6670 


0.7445 


0.9806 


0.9972 


0.9995 


1.0000 


0.1 


0.5517 


0.6687 


0.7462 


0.8683 


0.9421 


0.9820 


0.9982 


1.0003 


1.0009 


0.2 


0.5620 


0.6788 


0.7560 


0.8775 


0.9502 


0.9885 


1.0033 


1.0056 


1.0070 


0.5 


0.6593 


0.7714 


0.8402 


0.9571 


1.0230 


1.0615 


1.0849 


1.0910 


1.0939 


1.0 


0.9754 


1.0SS5 


1.1563 


1.2796 


1.3709 


1.4399 


1.4864 


1.4988 


1.5045 


2.0 


1.8362 


1.9938 


2.1113 


2.3336 


2.5129 


2.6508 


2.7439 


2.7685 


2.7799 


5.0 


4.5815 


4.9687 


5.2605 


5.8158 


6.2638 


6.6088 


6.8413 


6.9030 


6.9315 


10.0 


9.1629 


9.9373 


10.5209 


11.6315 


12.5276 


13.2176 


13.8060 


13.8060 


13.8629 



Fig. 9 — Semi-infinite slice with probe perpendicular to edge (or quarter- 
infinite slice with probe on a nonconducting edge). 
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results from Fig. 6 can be incorporated into the calculations, for we may 
let the correction divisor 

CD. = CD. (for nonconducting infinite slices) ± AC.D., (17) 

where the plus sign is used if the edge is nonconducting and the minus 
sign if the edge is conducting (the faces of the slice are assumed to be 
nonconducting in either case). Then 

1 1 



AC.D. = 



or 



where 



vV + 


M Vv 2 + i 

+ m [;1/(mVV + 


1)- 




■ m [in 


M fcvV + H)\ 


3ttV + *(mv^+ 

T ~r 74 


H) 


- N{nVv 2 + i) 




u = — and 

w 


v = 


L 

s 



(18) 



(19) 



(20) 

As shown in Fig. 8, L is the distance from the edge to the probe, s is 
the spacing of the probes, and w is the slice thickness. A table of values 
of the correction divisor is given in Fig. 8. 

The situation shown in Fig. 9 is similar to that of Fig. 8 in that the 
linear four-point probe is near the edge of a semi-infinite slice. But in 
this case the probe is perpendicular to the edge and L is the distance from 
the edge to the nearest point of the probe. Here again the correction 
divisor can be written in terms of the result for infinite slices, equation 
(17). It is found that 

ACD. = ± (- + I - - - I + M(y) + M(8) - M(a) - M(0)) (21) 

or 

ACD. = ~ (2I11 2| + N(a) + N(f3) - N(y) - N(8)) (22) 

2w \ a/3 / 



where 



a = (L + Ms)/w 

ft = (L + *As)/w 

(23) 

7 = (L + s)/w 

8 = (L + 2s)/w 
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Correction divisors can be compounded still further. For example, 
the correction divisor for the quarter-infinite slice with probe on the 
diagonal, shown in Fig. 10, may be written to make use of the results 
for Fig. 9. That is, 

CD. = CD. (Semi-infinite slice with probe perpendicular to (24) 
non-conducting edge, see Fig. 9) ±AC.D. 

Again, the plus sign is taken if both edges are nonconducting, the minus 
sign if both edges are conducting. If a, /3, 7, and 5 are redefined by 

1 



(25) 



a = — s/U + OS + L) 2 

fi = 2I V^+2s) 2 + (L + 3s) 2 

7 = ^ViXTJlTW 

8 = i- V(L + 3s)* + (L + S)> 

then AC.D. can be calculated by inserting equation (25) into equation 
(21) or (22) and doubling the result. One numerical example will be 




Fig. 10 — Linear probe on quarter-infinite slice. 



-J 
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given. If s = L = w, CD. = 2.1097 for nonconducting edges and 1.1665 
for conducting edges. 

It is evident that tabulation of correction divisors for situations with 
two or more degrees of freedom is rather tedious. Yet the computation 
of a correction divisor for a given set of dimensions may be very simple, 
requiring nothing more than the reading of several values of M or N 
from Fig. 3 or Fig. 4. 

In a four-point resistivity measurement it is not necessary that the 
points be equally spaced nor that they be colinear. Whatever the chosen 
arrangement, the correction divisor can be calculated from the M and N 
functions for any position of the probe on infinite, semi-infinite, or 
quarter-infinite slices. Numerous less frequently encountered situations 
may also be treated. 

The square arrangement of probe points illustrated in Fig. 11 is useful. 
When used on the surface of a semi-infinite solid, the resistivity is 
given by 

2tt6 AV 

P = P0 = 2~^V1 T <&> 




b/w 


0.1 


0.2 


0.5 


1.0 


2.0 


5.0 


10.0 


CD. 


1.0005 


1.004 


1.057 


1.344 


2.378 


5.916 


11.832 



Fig. 11 — Square arrangement of points — infinite slice. 
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where b is the length of the side of the square. When the square probe 
is used on infinite slices with nonconducting faces, the resistivity cal- 
culated from equation (26) must be corrected by dividing by 



CD. = 1 + 



or 



(2 - V2)io 

b 



[-( 



h 



- M I — 

,2w 



C - D - " (2^W [ M + N (k) ~ N ^k)] 



(27) 



(28) 



Some values of this correction divisor are given in Fig. 11. For thin slices, 



2twAV nn „ AV 
_. n = 9.06w -^r 
71n2 I 



(w«lb) 



(29) 



BASE RESISTANCE OF A POINT CONTACT TRANSISTOR AND TWO-POINT 

RESISTIVITY MEASUREMENTS 

The base resistance r» of a transistor is defined by 

n = (dVJdh) It (30) 

where V t is the emitter voltage, I c is the collector current, and I t is 
the emitter current. The main contribution to r b for a wide-spaced 
point-contact transistor seems to be the ohmic resistivity of the semi- 
conductor. In close-spaced point-contact transistors, the widening of 



COLLECTOR 



EMITTER 




-MAKE CONTACT 



s/w 


0.1 


0.2 


0.5 


1.0 


2.0 


C.F. 


0.931 


0.862 


0.667 


0.401 


0.118 



Fig. 12 — Base resistance of a point-contact transistor. 
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the collector space-charge region and the modulation of conductivity by 
carrier injection have important effects on / - 6 ; it is necessary to cal- 
culate the ohmic base resistance to evaluate these effects. The ohmic 
base resistance for a transistor made on an infinite slice with a base 
contact on the underside, as shown in Fig. 12, can be calculated from 
the M or X functions. It is found that 



!>- JL(l-Iln2+A \m(±)-m(± 

2ts \ w 2w\ \2wJ \4w 



(30) 



Hi 

[»(£)- »(b)] 

These results can be expressed in terms of a correction factor: 



* = /- 



n = ^- X C.F. (32) 

where p/2irs is the ohmic base resistance for a transistor constructed on 
a semi-infinite solid with base contact at infinity. Some values of the 
correction factor are given in Fig. 12. 

In measuring the resistivity of a slice which has been provided with a 
conducting coating on one face, it is better to use two points in the 
geometry shown in Fig. 12 than to use a four-point probe, for the me- 
chanical reason mentioned previously. The two-point probe also gives 
better localization of the measurement and reversal of the current 
affords a valuable check on the assumption of ohmic flow. Obviously, 
the correction divisors for correcting such a resistivity measurement 
for slice thickness are equal to the base resistance correction factors 
given in Fig. 12. 

FURTHER APPLICATIONS OF THE TABULATED FUNCTIONS 

In the previous examples, the potential was required only along a 
line perpendicular to a line of point sources and passing through one of 
the sources. The tabulations are of more general usefulness. For example, 
[Af(A) — 2M(2\)]q/a gives the potential of a line of point sources along 
a perpendicular bisector of the line between adjacent sources (this po- 
tential is zero when X = 0; the potential at infinity is — *=.). 

The potential of a plane grid of equal point sources is easily obtained 
from the tabulations. Suppose that the sources have a regular spacing a 
in one direction and b in the perpendicular direction. There is no loss of 
generality in assuming that b is greater than or equal to a. The potential 
at a point a distance d from the plane of sources and on a perpendicular 
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to the plane erected at one of the sources is 

V = I RM (34) 

a 

where X = d/a, k = b/a, and 

AW-i-JfW-Jjfg) 



(35) 



- 2 £ [M(VX 2 + nW) - M(nk)] 



Equation (4) enables one to transform equation (35) into the following 
rapidly converging expressions: 

sinh irX/k 

(36) 



R M = - - MOO - 2 In 



- 2 £ [#(»*0 - iV(Vx 2 + ^ 2 )] 



B ft (X) = tf(\) - 2 In 



sinh ?rX//b 
vk/k 



(37) 



- 2C - 2 E [#(nfc) - iV(VX 2 + n 2 fc 2 )] 

R-l 




Fig. 13 — Linear probe on square filament. 
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Equation (36) is best for small values of X, equation (37) for large values. 
In either case, one need consider only n = 1 to get five-place accuracy. 
The potential of a plane grid of sources may he used to calculate the 
correction divisor for four-point resistivity measurements on filaments. 
The case of square filaments is illustrated in Fig. 13; the correction divisor 
for this case is 



CD. = * [* (*) - R, (■ 

W |_ \10/ \1 



(38) 



For example, when s/w = %, CD. = 3.10. This value agrees very well 
with some experiments in which the reading of a four-point probe on a 
long germanium filament was compared with the resistivity determined 
from the potential gradient in the filament when current flowed through 
its entire length. 

The usefulness of the tabulations for three-dimensional arrays of point 
sources was tested by calculating the Madelung constant a for the 
Coulomb energy of an ionic crystal of the sodium chloride type. The 
calculation involves summing the potentials of a number of lines of 
sources of alternating sign; these potentials are obtained in the same 
way as for the calculation of the base resistance of a point-contact 
transistor. Nine terms must be considered to make full use of the ac- 
curacy of the tables. It is found that 

a = 2 In 2 + 4 \n Q - N fa@\ - 2N{\) + 2N fe^ 



- N it) - 2N (^) + 2N (m (39) 



-# (| V2) -4iV(V5)l 



and a value 1.74755 is obtained. The correct value is 1.747558. 3 

Plane boundaries between media of different but finite resistivities 
(or different dielectric constants or different thermal conductivities) can 
be treated by the method of images. When infinite image systems arise 
in this sort of problem, the tabulations given in this article can be 
applied. 

CONCLUSION 

The potential of a line of point sources can be obtained from either 
of the functions M and .V, which have been defined and tabulated. In 
addition to the tabulated functions, only elementary functions are re- 
quired for this purpose. 
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The functions M and N represent the potentials of definite source 
systems. This fact enables one to visualize the construction of solutions 
to potential problems involving lines of point sources. 
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APPENDIX 

From Fig. 1 , the superposition principle, and the usual q/r potential 
of a point source, it may be seen that 



*« - 2 £ 6 - vm) 



(40) 



According to the binomial theorem, 

1 _i . / 1\ -;u2 . ( \\( 3\ 1 „-5, 



= n 



+ - «-V + - - U.-V... (41) 



V^+T 2 'V 2/ ' \ 2 A 2/ 2! 

so that 



M(\) = -2 £ 



B>-v+HX-IH"^--> 



or 



M(X) = f(3)X 2 - Ht (5)X 4 + • • • , (43) 

where 

M = E n 9 (44) 



7. = 1 



is the Riemann zeta function. Equation (43) indicates a convenient 
approximation to M(\) for very small X: 

M(\) tt f(3)X 2 = 1.202X 2 (45) 

The binomial expansion, equation (41), converges more rapidly, the 
higher the value of n. It is, therefore, worthwhile to split up the sum in 
equation (40) into two parts. The first part will be a sum of the terms 
from n= 1 to n = m, calculated with the original expression for the 
summand. The second part will be a summation from » = m + 1 to 
n = co of the binomial expansion of the summand. It is expedient to 



THE POTENTIALS OF INFINITE SYSTEMS OF SOURCES 



127 



take in = 1 for X < 1 . Thus, 

M{\) = 2 - 2(1 + \ 2 r 1/2 + A 3 \ 2 - %A h \ A + %Aik* 

-35^ gX « + Gy 128 A u \ 10 - 23 Mi2^i3X 12 + 
where 



(46) 



a, = m - i = Z 



(47) 



Table II gives numerical values of the coefficients of equation (46). 5 
Table I is based mainly on calculations with in = 3. These were done 
at the Laboratories' IBM installation. 

Table II — Numerical Values of Coefficients in Equation (46) 

M{\) = 2 - 2(1 + X 2 ) -1 ' 2 + .20205 69032X 2 - .02769 58163X 4 

+ .00521 82984X 9 - .00109 83398X 8 + .00024 32335X 10 
- .00005 53648X 12 + ••• 



Table III — Some Values 


of the Derivative Af'(X) 


N 


-l/'(X) 


X 


-1/'(X> 


0.05 


0.11982 


1.30 


0.95001 


0.10 


0.23734 


1.35 


0.93512 


0.15 


0.35040 


1.40 


0.92004 


0.20 


0.45709 


1.45 


0.90489 


0.25 


0.55586 


1.50 


0.88975 


0.30 


0.64556 


1.55 


0.87471 


0.35 


0.72545 


1.60 


0.85982 


0.40 


0.79519 


1.65 


0.84513 


0.45 


0.86714 


1.70 


0.83068 


0.50 


0.90467 


1.75 


0.81649 


0.55 


0.94528 


1.80 


0.S0259 


0.00 


0.97736 


1.85 


0.78900 


0.G5 


1.00174 


1.90 


0.77568 


0.70 


1.01926 


1.95 


0.76270 


0.75 


1.03077 


2.00 


0.75003 


0.80 


1.03709 


2.05 


0.73768 


0.85 


1.039(H) 


2.10 


0.72564 


0.90 


1.03719 


2.15 


0.71391 


0.95 


1.03229 


2.20 


0.70249 


1.00 


1.02487 


2.25 


0.69136 


1.05 


1.01541 


2.30 


0.68053 


1.10 


1.00431 


2.35 


0.66999 


1.15 


0.99195 


2.40 


0.65972 


1.20 


0.97862 


2.45 


0.64973 


1.25 


0.96457 


2.50 


0.64000 
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Table IV — Some Values of the Derivative JV'(X) 



X 


-A'(X) 


X 


-A"(X) 


1.00 


0.02487 


1.55 


0.00063 


1 05 


0.01766 


1.60 


0.00046 


1.10 


0.01257 


1.65 


0.00033 


1.15 


0.00896 


1.70 


0.00024 


1.20 


0.00640 


1.75 


0.00019 


1.25 


0.00457 


1.S0 


0.00012 


1.30 


0.00326 


1.85 


0.00009 


1.35 


0.00233 


1.90 


0.00006 


1.40 


0.00167 


1.95 


0.00004 


1.45 


0.00120 


2.00 


0.00003 


1.50 


0.00086 







Although equation (46) certainly converges for all values of X less 
than 2, it is not efficient for X > 1. For X > 1, it is better to work with 
the function N. According to Madelung, 



N(\) = 2t £ iH w (i2im\) 



(48) 



where iH ll) (ix) is the Hankel function tabulated in Jahnke & Emde. 
The functions M(X) and N(\) are tabulated in Table 1. M(K) and 
iV(X) are connected by equation (4). Because of the overlapping regions 
of X-values for which the expressions for M and N converge with reason- 
able rapidity, it is possible to tabulate both M and N for all values of 
X. However, for satisfactory graphical interpolation, it is necessary to 
use the smaller of the two functions, so that M is given in Fig. 3 for small 
values of X and N is given in Fig. 4 for large values of X. Table III gives 
some values of the derivative, Af'(X). The field E of a line of point charges 
can be calculated from M' in accordance with the equation obtained by 
differentiating equation (3) with respect to d. Thus 



E 



.ifi 

a 2 |_X 2 



+ M'(X) 



(49) 



The derivative N'{\) is tabulated in Table IV. 
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